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ABSTRACT 


'‘Teledyne  Brown  Engineering  (TBE)  has  developed  high-fidelity 
exoatmospheric  time-varying  radar  cross  section  (RCS)  models  for  use 
in  the  Warning  Information  Correlation  (WIC)  threat  development  effort. 
These  models  have  been  developed  for  the  objects  of  six  USSR  ICBM  and 
SLBM  weapon  systems.  This  report  describes  the  five  TBE  developed  time- 
varying  RCS  modeling  algorithms  that  were  utilized  for  these  models. 

Part  of  the  validation  of  these  models  consisted  of  comparisons  of  selected 
discrimination  signatures  calculated  from  model  and  recorded  data.  The 
definition  of  these  signatures  is  also  presented./ - 
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1.  INTRODUCTION 
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I 


TBE  was  tasked  under  Contract  DASG60-79-C-0094  to  develop  high- 
fidelity  ex oatmo spheric  time-varying  RCS  models  of  the  objects  of  six 
USSR  ICBM  and  SLBM  systems  at  L-band  and  UHF  frequencies.  These  models 
have  been  developed  and  validated  through  analysis  and  comparisons  with 
recorded  COBRA  DANE  and  ARIS  RCS  data  on  USSR  missile  system  tests.  These 
data  are  maintained  in  rapid-access  computer  files.  Currently  the  L-band 
RCS  data  base  contains  5629  tracks  from  248  COBRA  DANE  missions  covering 
21  different  missile  systems.  These  L-band  data  represent  a  total  of  over 
336  hours  of  RCS  observations.  The  UHF  data  base  currently  includes  473 
histories  representing  approximately  7  hours  of  ARIS  RCS  observations 
covering  15  different  missile  systems. 

Five  algorithms  were  developed  by  TBE  to  simulate  the  sets  of  observed 
RCS  data  characteristics  of  objects  to  be  modeled.  These  algorithms  utilize 
discrete  inputs  from  a  tunable  variable  space  defined  for  each  modeled 
object.  These  tunable  parameters  were  generally  specified  in  terms  of 
probability  distributions  for  each  model  developed  to  provide  for  observed 
mission-to-mission  variations  in  the  recorded  RCS  data.  Section  2  of  this 
report  presents  a  brief  description  of  the  time-varying  RCS  modeling  algorithms 
that  were  developed.  A  detailed  description  of  each  algorithm  Is  given 
in  the  appendices.  An  identification  of  the  algorithm  utilized  for  each 
object  of  the  USSR  weapon  systems  that  have  been  modeled  is  given  in 
Reference  1. 

Section  3  provides  a  discussion  of  the  procedures  used  for  RCS 
model  validation.  Part  of  this  procedure  consisted  of  comparing  the 
values  of  discrimination  signatures  calculated  from  the  recorded  and  the 
model  data.  A  definition  of  these  signatures  is  provided.  The  potential 
effects  of  random  noise  on  these  signatures  is  also  discussed. 


2.  RCS  MODEL.  ALGORITHMS  DEVELOPMENT 

During  the  process  of  developing  time-varying  RCS  models  of  the 
various  objects  of  six  USSR  ICBM  and  SLBM  weapon  systems,  various  algorithms 
were  developed  to  simulate  the  observed  RCS  time  structure  and  RCS  discrimina¬ 
tion  signatures  of  the  recorded  data  at  L-Band  and  UHF  frequencies.  The 
RCS  modeling  proceeded  concurrently  with  the  TBE  development  of  flight 
models  for  the  same  weapon  systems.  As  objects  were  identified  by  the 
flight  modeling  efforts,  the  associated  RCS  data  were  analyzed.  Differing 
RCS  data  structures  observed  for  the  different  object  types  has  led  to 
the  development  of  several  algorithms. 

These  algorithms  were  separately  documented  as  they  were  developed 
by  TBE.  These  documents  have  been  edited  and  revised  to  reflect  the  current 
form  of  the  algorithm  and  are  included  as  Appendices  A  through  E.  A  brief  descrip 
tion  of  each  algorithm  is  given  in  the  following  sections. 

2.1  SPB  ALGORITHM 

The  Static  Pattern  Based  (SPB)  algorithm  has  been  employed  to 
develop  RCS  models  for  large,  well  defined,  roll-symmetric  objects  (e.g. 
cylinders,  cones,  etc)  with  well  defined  body  motion.  This  algorithm 
provides  viewing  and  body  motion  sensitivities.  Observed  RCS  data  for 
these  objects  are  largely  cyclically  repetitive  as  a  function  of  time. 

An  RCS  static  pattern  (RCS  in  dBsm  vs  aspect  angle)  is  used  in  this  algorithm 
to  provide  an  RCS  estimate  as  a  function  of  the  aspect  angle  which  is 
computed  from  the  modeled  body  motion  and  the  radar  viewing  parameters  at 
each  model  simulation  time.  The  RCS  static  patterns  for  this  algorithm 
are  based  upon  analytic  calculations  from  intelligence  estimates  of  the 
object  size  and  shape  as  well  as  the  observed  RCS  data.  The  SPB  algorithm 
is  documented  in  Appendix  A. 
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2.2 


MSPB  ALGORITHM 


The  Modified  Static  Pattern  Based  (MSPB)  algorithm  is  similar 
to  the  SPB  algorithm,  but  also  models  spin  motion  for  non-roll  symmetric 
objects.  This  algorithm  also  supports  the  modeling  of  scattering  from 
objects  composed  of  multiple  bodies.  The  MSPB  algorithm  models  both  the 
main  body  motion  and  the  relative  motion  of  any  additional  scatterers. 

The  model  parameters  include  motion  parameters  which  are  trajectory  de¬ 
pendent  as  well  as  motion  parameters  which  are  trajectory  independent. 

The  MSPB  algorithm  is  detailed  in  Appendix  B. 

2.3  CORS  ALGORITHM 

TBE  developed  the  Correlated  Statistical  (CORS)  algorithm  for 
objects  whose  RCS  histories  are  essentially  random  in  nature.  This  model 
does  not  include  any  body  motion  or  radar  viewing  parameters.  The  CORS 
algorithm  simulates  the  observed  RCS  data  by  computing  correlated  random 
draws  from  a  constrained  Gaussian  distribution  in  dBsm  space.  The  correla¬ 
tion  function  is  modeled  as  an  exponential  time  decay.  The  mathematical 
description  of  the  CORS  algorithm  is  provided  in  Appendix  C. 

2.4  SARP  ALGORITHM 

The  Sinusoids  and  Random  Process  (SARP)  algorithm  was  developed 
to  model  RCS  characteristics  that  contain  significant  time  structure  that 
cannot  be  associated  with  definable  object  motion  characteristics.  Spectral 
analysis  of  these  RCS  data  show  a  pronounced  contribution  from  a  limited 
sequence  of  sinusoids  in  addition  to  a  broad-band  background.  The  SARP 
model  simulates  this  observed  structure  by  summing  the  contribution  from 
a  specified  number  of  harmonically  related  sinusoids  (in  dBsm  space)  and 
a  constrained  uncorrelated  Gaussian  background  process.  The  number  of 
sinusoids,  their  amplitudes  and  phases  are  defined  from  the  power  spectra 
of  the  RCS  data.  The  SARP  algorithm  and  the  required  model  parameters  are 
described  in  Appendix  D. 
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2.5 


SARP2  ALGORITHM 


The  SARP2  algorithm  was  developed  for  one  object  whose  spectral 
analysis  showed  only  a  limited  number  of  non-harmonic  frequency  components. 
Functionally,  the  SARP2  algorithm  is  nearly  identical  to  the  SARP  algorithm 
(Appendix  D) .  The  primary  difference  between  SARP  and  SARP2  is  that  SARP2 
uses  explicitly  defined  sinusoids  that  are  not  constrained  to  be  harmonic. 
Thus  the  SARP2  algorithm  requires  the  frequency,  amplitude  and  phase  to 
be  defined  for  each  sinusoid. 

2.6  MSARP  ALGORITHM 

The  Modified  Sinusoids  and  Random  Process  (MSARP)  algorithm  was  de¬ 
fined  for  objects  whose  RCS  spectral  data  possessed  many  non-harmonically 
related  sinusoids.  These  objects  are  more  irregular  than  those  modeled 
by  SARP.  Since  the  MSARP  algorithm  does  not  explicitly  include  motion 
parameters,  this  algorithm  also  does  not  support  radar  viewing  dependencies. 
The  MSARP  algorithm  is  described  in  Appendix  E. 
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3,  RCS  MODEL  VALIDATION 

The  time-varying  RCS  model  for  an  object  consists  of  the  specifica¬ 
tion  of  the  time-varying  algorithm  and  the  required  probability  distribu¬ 
tions  for  the  algorithm  input  parameters.  Validation  of  a  mode]  was  comprised 
of  two  aspects.  The  initial  aspect  was  accomplished  by  visual  comparisons 
of  the  structure  of  the  real  and  simulated  time  histories.  Recorded  RCS 
time  histories  and  those  produced  from  the  model  should  show  similar  structure 
and  basic  characteristics.  A  correspondence  should  be  noted  in  the  troughs 
and  peaks  with  the  same  approximate  periodicity.  Rather  than  a  one-to-one 
correspondence  between  an  RCS  history  and  a  model  realization,  the  model 
realizations  should  be  representatives  of  the  basic  characteristics  observed 
in  the  data  of  a  particular  object.  The  second  aspect  of  model  validation 
was  a  comparison  between  selected  RCS  discrimination  signatures  calculated 
from  both  the  recorded  data  and  the  time-varying  RCS  model.  The  RCS 
signatures  selected  for  model  validation  are  identified  in  Section  3.1. 

The  effects  of  random  noise  on  these  signatures  are  given  in  Section  3.2. 

3.1  RCS  DISCRIMINATION  SIGNATURES 

The  selected  RCS  signatures  include  the  RCS  signatures  currently 
employed  by  the  PARCS  for  characteristic  object  processing,  the  signatures 
used  by  the  PAVE  PAWS  automatic  object  classification  subfunction  and  the 
signatures  used  by  COBRA  DANE  (except  the  coefficient  of  variation).  The 
power  discrimination  ratio  which  had  been  proposed  for  EPARCS  was  also 
selected  as  well  as  two  RCS  signatures  previously  proposed  by  TBE. 

The  calculated  RCS  discrimination  signatures  are: 

•  peak  RCS  (PRCS) 

•  average  RCS  (ARCS) 

•  maximum  absolute  pulse  to  pulse  RCS  difference  (MPD) 

•  average  absolute  pulse  to  pulse  difference  (APD) 

•  power  discrimination  ratio  (PDR) 

•  RCS  variance  (VRCS) 
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•  scintillation  index  (SCIN) 

•  local  variance  (LVS) 

•  local  residual  (LRS) 

All  signatures  except  the  PDR  are  calculated  in  dBsm  space.  The  PDR  is 
calculated  in  M2  space.  The  PDR  is  essentially  the  ratio  of  low  frequency 
content  to  the  total  RCS  observation.  The  SCIN  can  be  considered  as  the 
ratio  of  the  high  frequency  to  low  frequency  content.  The  LVS  is  the 
average  variance  calculated  over  sets  of  four  RCS  observations.  The  LRS 
is  similar  to  LVS  except  the  mean  square  distance  is  calculated  from  a 
least  squares  straight  line  fit  to  the  four  points.  The  MPD,  APD,  LVS, 
and  LRS  are  calculated  only  when  the  time  between  RCS  observations  is 
less  than  0.3  seconds.  The  LVS  and  LRS  require  at  least  4  consecutive 
observations  meeting  this  condition.  The  PDR  requires  data  spaced  at  a 
maximum  of  1.1  seconds.  These  selected  RCS  discrimination  signatures,  a 
brief  summary  of  calculation  procedures  and  restrictions  are  summarized 
by  Table  3-1. 

3.2  NOISE  EFFECTS  ON  DISCRIMINATION  SIGNATURES 

Recognizing  the  potential  for  the  effects  of  noise  on  the  signatures 
calculated  from  the  recorded  RCS  data,  an  analysis  was  performed  to  illustrate 
the  magnitude  of  these  effects.  This  analysis  consisted  of  superimposing 
zero-mean  Gaussian  random  noise  constrained  to  ±3o  limits  on  various  constant 
ARCS  levels.  This  process  was  sampled  at  a  24  Hz  rate  for  a  period  of  50 
seconds  and  the  RCS  signatures  were  calculated  from  these  1200  samples. 

The  results  of  this  exercise  show  that  only  the  ARCS  signatures  were  un¬ 
affected  by  the  Gaussian  process.  Many  of  the  resulting  signatures  could 
be  approximated  from  the  standard  deviation  (o^)  of  the  simulated  Gaussian 
noise  by: 

PRCS  =  ARCS  +  3  o„ 

N 

MPD  =  4.8  * 

N 

APD  =  1.1  *  CL. 

N 

VRCS  =  aN2 

LVS  =  0.72  a.  2 
N 

LRS  =  0.49  a  2 
N 
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Corresponding  results  for  the  PDR  and  SC1N  signatures  are  shown  by  Figures 
3-1  and  3-2. 

This  exercise  indicates  that  any  zero-mean  Gaussian  noise  contained 
in  the  real  RCS  data  would  tend  to  increase  the  calculated  PRCS,  MPD,  APD, 
VRCS ,  LVS  and  LRS  signatures.  The  effects  on  SCIN  and  PDR  are  more  de¬ 
pendent  on  the  ARCS  and  the  ARCS  signature  itself  is  not  influenced  by 
zero-mean  noise.  These  effects  of  noise  on  the  signatures  were  considered 
when  tuning  an  RCS  model  to  obtain  signatures  that  correspond  to  those 
calculated  from  the  recorded  data. 
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*N  Pulses  assumed  with  N-l  pulse  differences  and  N/4  four  pulse  groups. 


TABLE  3-1.  RCS  DISCRIMINATION  SIGNATURES,  PROCEDURES,  AND  RESTRICTIONS  (concluded) 
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FIGURE  3-1.  EFFECTS  OF  RANDOM  NOISE  ON  PDR  WITH  CONSTANT  ARCS  LEVELS 
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ARCS  =  0 


cr|y  (dBsm) 


FIGURE  3-2.  EFFECTS  OF  RANDOM  NOISE  ON  SC  IN  WITH  CONSTANT 
ARCS  LEVELS 
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APPENDIX  A 


STATIC  PATTERN  BASED  (SPB)  TIME 
VARYING  RCS  ALGORITHM 


By 

W.  B.  Wellman 
D.  C.  Johnson 


1.  INTRODUCTION 


TBE  is  required  to  develop  exoatmospheric  time  varying  L-band 
and  UHF  RCS  models  for  objects  of  several  missile  systems  under  Contract 
DASG60-79-C-0094 .  A  static  pattern  based  (SPB)  time  varying  RCS  genera¬ 
tion  method  is  one  procedure  for  simulating  RCS  histories  that  have 
significant  time  structure.  TBE  has  developed  an  SPB  RCS  generation  model 
for  this  purpose. 

Section  2  presents  a  description  of  the  SPB  RCS  model.  The 
model  parameters  including  those  which  define  the  physical  motion 
characteristics  imposed  on  the  object  during  flight  are  given. 
Representative  model  output  for  several  motion  and  viewing  conditions 
is  presented. 
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2.  STATIC  PATTERN  BASED  RCS  GENERATION  MODEL 


2.1  INTRODUCTION 

The  RCS  scattering  characteristics  of  objects  that  are  roll 
symmetric,  e.g.  cylinders,  cones,  etc.,  can  be  represented  by  a  single 
RCS  static  pattern  for  a  specified  combination  of  radar  transmit  and 
receive  polarizations.  The  static  pattern  is  the  RCS  amplitude  as  a 
function  of  the  viewing  angle  relative  to  the  body  centerline.  This 
angle  is  called  the  aspect  angle.  If  the  position  and  orientation  of 
the  object  relative  to  the  radar  at  a  given  transmission  time  can  be 
determined,  the  aspect  angle  at  which  the  radar  is  viewing  the  object 
can  be  calculated.  Then  the  RCS  amplitude  from  the  static  pattern 
at  the  calculated  aspect  angle  is  the  RCS  for  that  transmit  time. 

Modeling  the  body  motion  as  the  object  flies  along  its  trajectory  and 
computing  the  aspect  angle  and  then  the  RCS  for  each  transmission  time 
produces  a  time  varying  RCS  for  an  object  as  viewed  from  a  specified 
radar  location. 

Thus  for  this  technique,  a  model  of  the  object  motion  must  be 
developed  and  the  equations  for  computing  the  aspect  angle  to  the  object 
from  a  given  radar  location  must  be  derived.  The  following  sections 
define  the  object  motion  model  and  viewing  geometry  relationships 
between  the  radar  and  object  orientation.  Then  in  Section  2.5  the  model 
implementation  is  described  and  example  time  varying  RCS  histories  are 
given  in  Section  2.6. 

2.2  BODY  MOTION  MODEL 

The  body  motion  is  defined  relative  to  a  body  motion  axis.  The 
body  motion  axis  is  defined  as  a  reference  axis  extending  through  the 
body's  center  of  gravity  at  the  position  of  the  body  on  its  ballistic 
trajectory  path.  The  orientation  is  defined  by  an  angle  t p  relative  to 
a  perpendicular  to  the  trajectory  plane  and  by  an  angle  T  from  a  reference 
axis  in  the  trajectory  plane  as  shown  in  Figure  1(a) . 


r  .  —  ■  1 

I 

i  s 

In  this  model  two  types  of  body  motion  are  included.  First,  i 

there  is  precession  motion,  which  is  the  fundamental  cyclical  motion  of 
the  body  about  the  motion  axis.  This  motion  is  assumed  to  occur  with 
the  body  centerline  moving  at  an  angular  rate  A  on  the  surface  of  a 
cone  P  degrees  from  the  motion  axis  as  shown  in  Figure  1(b).  An  object 
purely  tumbling  in  its  trajectory  plane  would  correspond  to  ip  =  0  and 
P  =  90°. 

The  second  type  of  body  motion  modeled  is  nutation,  which  is  an 

additional  precession  motion  of  the  body  about  the  mean  body  centerline 

axis.  This  motion  is  defined  by  a  angle  N  from  the  mean  body  centerline 

and  an  angular  rate  D  as  shown  in  Figure  1(c).  There  are  two  constraints 

on  the  angles  N.  First  N  is  restricted  to  be  less  than  P.  Second  there 

is  a  joint  constraint  between  the  allowed  value  of  N  and  the  values  of 
•  •  *  * 

D  and  A.  N  is  constrained  to  be  less  than  the  maximum  of  [15.6845  -  D/A]  x 

0.6588  and  [51.3543  -  D/A]  x  0.0627.  This  latter  constraint  which  was 

empirically  derived,  is  necessary  to  prevent  unrealistic  combination  model 

j  of  these  parameters. 

I:  2.3  RADAR  VIEWING  GEOMETRY 

[  For  the  SPB  RCS  model,  the  change  in  the  radar  line  of  sight  as 

,  the  object  moves  along  its  trajectory  path  is  modeled  in  an  approximate 

t 

1  manner.  At  an  initial  time  tQ,  e.g.,  when  or  before  the  initial  time  of 

the  desired  RCS  interval,  the  angle  0  between  the  radar  line  of  sight 
~  and  the  perpendicular  to  the  trajectory  plane  is  ietermined.  As  shown 

in  Figure  2,  the  radar  line  of  sight  vector  is  projected  into  the 
I  trajectory  plane.  For  convenience,  this  is  used  as  the  reference  axis 

for  the  angle  T  defined  in  Figure  1(a) .  The  angle  (j)  between  the  j 

radar  line  of  sight  and  the  body  motion  axis,  which  will  be  required  . 

1  in  Section  2.4,  is  also  identified  in  Figure  2. 

t 

The  change  in  the  viewing  angle  as  the  object  moves  along  its  j 

trajectory  path  is  modeled  by  assuming  that  9  changes  at  a  constant  rate 
9.  Thus  9(t)  =  0O  +  0 (At) ,  where  80  is  the  viewing  angle  at  time  t0  and  j 

At  =  t  -tQ.  | 
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ASPECT  ANGLE  DETERMINATION 


The  aspect  angle  will  first  be  determined  relative  to  the  body 

centerline  motion,  i.e.,  the  precession,  and  then  this  equation  will  be 

modified  to  account  for  the  nutation.  From  Figure  3,  the  aspect  angle 

c  is  c(t)  »  cos  ^  (sin  (t)  sin  P  cos  (AAt)  +  cos  4>  (t)  cos  P) ,  where  $ 

is  the  angle  between  the  radar  line  of  sight  and  the  body  motion  axis  and 

A(At)  is  the  angular  amount  that  the  body  centerline  has  rotated  since 

t ime  t  . 
o 

from  Figure  4,  which  includes  the  nutation  motion,  the  angle 

from  the  body  motion  axis  to  the  body  centerline  is  the  angle  P  plus  the 

change  due  to  the  nutation.  This  approximate  change  in  P  due  to  the 

nutation  is  N  .  cos  (DAt) .  The  change  in  the  precession  rotation  is 
-1 

approximately  tan  (tan  N  sin  (DAt)  sin  P).  Thus,  including  the 
nutation,  the  equation  for  the  aspect  angle  becomes 

e(t)  =  cos-1  [sin  $(t)  sin  (P  +  N  .  cos  (DAt)) 

.  cos  (AAt  +  tan  ^(tan  N  .  sin  (DAt)/sin  P))  (1) 

+  cos  (J) ( t )  cos  (P  +  N  .  cos  (DAt))]  . 

To  complete  the  computation  of  e(t),  <J>(t)  is  required.  From 

Figure  2, 

<p(t)  =  cos  ^  [sin  ip  sin(9Q  +  9At)  sin  T 

+  cos  ip  cos  (90  +  9At))].  , 

2.5  MISSION-TO-MISSION  VARIABILITY  PARAMETERS 

Recorded  RCS  data  from  objects  for  which  the  SPB  algorithm  was 

used  have  shown  an  unexplained  variability  in  the  peak  RCS  from  mission 

to  mission.  To  model  this  variation  a  small  amplitude  slowly  varying 

sinusoid,  A„  sin  (2ttT/T  ),  is  added  to  the  computed  RCS  values.  These 
°  s 

two  parameters,  Ag  and  Tg,  given  in  Table  3,  complete  the  parameter  defin 
tlon  of  the  model.  The  variable  T  is  the  sum  of  the  relative  simulation 
time  and  an  initial  uniform  random  draw  from  the  interval  (0,  T  ). 


! 


2.6  MODEL  IMPLEMENTATION 

From  the  preceding  sections,  six  motion  and  two  mission-to-mission 
variability  parameters  for  the  model  have  been  defined.  These  are  given 
in  Tables  1  and  3.  The  use  of  this  model  requires  a  definition  of  each 
parameter.  In  general,  these  definitions  will  be  in  the  form  of  proba¬ 
bility  density  functions  to  account  for  the  variability  in  object  motion 
that  is  typically  observed  on  various  test  flights  of  the  same  object. 

To  complete  the  SPB  RCS  model,  the  static  pattern  is  required. 

This  can  be  either  an  analytical  expression  or  a  table  of  RCS  amplitude 
as  a  function  of  the  aspect  angle.  In  the  latter  case,  the  RCS  is 
obtained  by  linear  interpolation. 

The  use  of  this  model  requires  the  definition  of  the  radar 
viewing  condition  parameters  given  in  Table  2.  The  value  of  9q  can  be 
approximated  or  calculated  from  the  relationship  between  the  radar  and 
the  object  on  its  ballistic  trajectory  at  time  t  . 

After  the  model  and  radar  viewing  parameters  have  been  established 

for  t  ,  the  use  of  the  equations  (1)  and  (2)  with  At  =  t  -t  provides  the 
o  o 

value  of  e  at  time  t.  The  RCS  at  time  t  is  determined  for  the  calculated 
value  of  t. 

2.7  SPB  RCS  MODEL  EXAMPLES 

Figure  5  shows  an  analytical  static  pattern  of  RCS  amplitude  in 

dBsm  as  a  function  of  aspect  angle  in  degrees.  The  pattern  was  generated 

2 

from  RCS  =  (sin  x/x)  ,  where  x  =  10.1664  cos  ft. 

Figure  6  and  7  show  examples  of  the  SPB  RCS  model  output  using 
the  static  RCS  pattern  in  Figure  5  and  with  a  10-Hz  sample  rate.  Figure  6 
shows  the  effect  of  changes  in  the  radar  location  and  precession  motion 
on  the  time  varying  RCS  history.  Figure  7  shows  the  effects  of  nutation 
motion  with  the  radar  in  the  trajectory  plane  and  the  object  tumbling  in 
the  plane  of  the  trajectory.  The  effects  of  precession  motion  and  radar 


location  changes  result  in  varying  the  region  of  the  static  RCS  pattern 

within  which  sampling  takes  palce.  The  addition  of  nutation  motion 

results  in  increased  scintillation  in  the  time  varying  RCS.  In  Figure  6 

and  7,  A  and  0  set  to  zero  because  the  time  intervals  were  too  short  to 
s 

notice  an  appreciable  effect  on  the  time  varying  RCS  due  to  trajectory 
and  radar  viewing  geometry  changing  with  time. 
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TABLE  1.  SPB  RCS  MODEL  MOTION  PARAMETERS 


SYMBOL 

DESCRIPTION 

Angle  between  body  motion 
axis  and  perpendicular  to 
trajectory  plane  (Deg) 

T 

Angle  in  trajectory  plane 
between  projection  of  body 
motion  axis  and  radar  line 
of  sight  at  time  t0  (Deg) 

p 

Precession  Angle  (Deg) 

T  _  360 

p  A 

Precession  Period  (sec) 

N 

Nutation  Angle  (Deg) 

D 

Nutation  Angular  Rate 
(Deg/sec) 

TABLE  2.  RADAR  VIEWING  CONDITION  INPUT  FOR  SPB 

RCS  MODEL 


SYMBOL 

DESCRIPTION 

^0 

Angle  between  radar  line 
of  sight  and  perpendicular 
to  trajectory  plane  at 
time  t0  (Deg) 

9 

Angular  rate  of  change 
of  9  (Deg/sec) 
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TABLE  3.  SPB  RCS  MODEL  MISSION-TO-MISSION 
VARIABILITY  PARAMETERS 


SYMBOL 

DESCRIPTION 

A 

s 

Sinuisoid  Amplitude  (dBsm) 

T 

s 

Sinuisoid  Period  (sec) 

.-9 


FIGURE  3.  BODY  PRECESSION  MOTION  GEOMETRY 


ASPECT  ANG 
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A  =  18  deg/sec 
N  =  0  deg 
D  =  0  deg/sec 
0  =  90  deg 


ip  =  1 0  deg 
t  =  180  deg 
P  =  90  deg 
A  =  18  deg/sec 
N  =  0  deg 
D  =  0  deg/sec 
0  =  50  deg 


ip  =  20  deg 
T  =  0  deg 
P  =  60  deg 
A  =  18  deg/sec 
N  =  0  deg 
D  =  0  deg/sec 
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1.  INTRODUCTION 


Teledyne  Brown  Engineering  (TBE)  is  required  to  develop  exo- 
atmospheric  time  varying  L-band  and  UHF  radar  cross  section  models  for 
objects  of  several  missile  systems  under  Contract  DASG60-79-C-0094.  One 
relevant  class  of  objects  requires  a  particular  spatial  orientation 
during  flight  and  is  not  in  general  roll-symmetric.  TBE  has  developed  a 
modified  static-pattem-based  (MSPB)  RCS  generation  method  for  modeling 
such  objects.  Section  2  of  this  report  describes  the  general  nature  of 
the  objects  being  considered  and  the  relative  motions  such  objects  are 
likely  to  possess.  The  object  and  object  motion  parameters  are  trajectory 
independent  variables.  The  spatial  orientation  of  the  object  and  the 
relationship  of  this  orientation  to  the  motion  of  the  object  is  detailed  in 
Section  3.  The  parameters  associated  with  the  particular  spatial  orienta¬ 
tion  required  by  the  object  during  flight  are  trajectory  dependent  variables. 
Section  4  outlines  the  method  used  to  tie  the  radar  location  to  the  spatial 
orientation  of  the  object.  The  parameters  necessary  to  affect  this  charac¬ 
terization  are  viewing  dependent  variables.  Section  5  then  presents  a 
summary  which  describes  the  method  by  which  all  the  parameters  are  brought 
together  to  form  a  time-varying  RCS  model.  Finally,  Section  6  presents 
an  example  implementation  of  the  model.  Much  of  the  mathematics  used  in 
this  model  is  detailed  in  Appendix  A,  and  references  to  equations  from  the 
appendix  are  prefixed  with  an  A. 
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2.  TRAJECTORY  INDEPENDENT  VARIABLES 


2.1  OBJECT  DEFINITION  (BODY  COORDINATE  SYSTEM) 

When  the  dimensions  of  the  object  in  terms  of  the  viewing  radar 
wavelength  are  large,  the  general  technique  for  estimating  the  object's 
RCS  consists  of  establishing  a  model  in  which  the  RCS  is  due  to  a  number 
of  component  scatterers.  Once  a  model  has  been  established  for  the 
various  component  scatterers,  the  relative  magnitudes  of  the  scattered 
fields  from  these  components  are  determined,  and  their  contributions  are 
then  added  with  the  proper  phase  to  provide  an  estimate  of  the  RCS  of  the 
entire  object.  Such  a  procedure  must  be  carried  out  for  all  possible, 
non-redundant  aspect  angles  to  obtain  a  complete  RCS  static  pattern  for 
the  object.  In  the  case  of  roll-symmetric  objects,  the  symmetry  permits 
tabulation  of  complete  RCS  information  in  a  relatively  short  listing. 

Such  symmetry  is  made  use  of  in  TBE's  static  pattern  based,  time  varying 
RCS  generation  model  described  in  Appendix  A.  The  fact  that  all  objects 
are  not  roll-symmetric  has  necessitated  the  development  of  the  MSPB  model. 

Figure  1  represents  an  asymmetric  object  with  a  number  of 
scattering  components  which  are  fixed  in  a  body  coordinate  system.  The 
origin  of  the  coordinate  system  is  located  at  the  center  of  gravity  of 
the  object,  and  the  scattering  component  most  closely  associated  with 
the  origin  is  termed  the  primary  scatterer.  The  vector  location  of  the 
ich  scattering  component  and  its  RCS  contribution  are  designated  S^  and 
Of,  respectively.  The  accuracy  of  the  MSPB  model  improves  as  the  magni¬ 
tude  of  each  S^  gets  larger  relative  to  the  radar  wavelength.  Also  note 
that  each  scattering  component  may  or  may  not  be  physically  connected  to 
other  scatterers  and  in  fact  may  move  about  in  some  random  or  determin- 
istic  fashion  in  the  body  coordinate  system.  That  is,  S^  may  be  a 
function  of  time. 

The  vector  which  points  to  the  viewing  radar  location  from  the 

A 

origin  of  the  body  coordinate  system  is  the  LOS  vector.  As  shown  in 
Figure  2,  this  vector  forms  a  polar  angle  0^  and  an  azimuth  angle  <J>^ 
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in  the  body  coordinate  system,  and  the  RCS  contribution  from  each 
scattering  component  is  a  function  of  either  one  or  both  of  these  angles. 
Namely , 


a 


i 


f  . 
1 


(1) 


The  form  of  the  function  depends  on  the  scattering  component  and  may 
range  from  a  static  pattern  look-up  table  or  a  simple  algorithm  to  some 
very  involved  analytic  expression  which  is  being  used  to  estimate  the  RCS 
of  the  complex  shape. 


The  proper  phase  of  each  scattering  component  (<3>i)  is  obtained  by 
referencing  the  phase  to  the  origin  of  the  coordinate  system  and  adding 
the  phase  delay  (A^)  of  each  scatterer  to  that  value.  That  is: 


<J>.  =  k  S.  •  LOS  +  A. 
11  l 


(2) 


where  k  =  2n/\  is  the  wave  number.  The  estimated  RCS  of  an  object  with 
M  scattering  components  is  given  by: 


RCS 


M  j  2<J> .  2 

is 7  e  1 

i-1  1 


(3) 


The  RCS  value  calculated  in  Equation  (3)  will  vary  with  time  because  of 
changing  values  of  0  and  <j>  •  These  angles  vary  because  the  object  both 

Li  Li 

flies  along  a  trajectory  and  moves  about  its  center  of  gravity. 

2.2  OBJECT  MOTION  (MOTION  AXIS  COORDINATE  SYSTEM) 

The  motion  of  a  rigid  body  in  the  absence  of  external  torques 
can  be  specified  by  three  independent  parameters  called  the  Eulerian 
angles  (a,  g,  y)  .  The  angular  momentum  vector  of  this  rigid  body  is 
constant  and  so  has  a  fixed  direction  in  space  which  in  this  report  is 
called  the  motion  axis  (MA) .  Spin  is  a  periodic  rotation  of  a  rigid 
body  about  its  own  centerline.  However,  as  the  body  spins,  an  observer 
fixed  in  the  body  coordinate  system  will  see  a  rotation  of  the  spin 


angular  velocity  vector  about  the  motion  axis.  This  second  rotation  is 
called  precession.  The  angle  between  the  body  centerline  and  the  motion 
axis  is  called  the  coning  angle  and  is  itself  subject  to  a  periodic 
change  called  nutation. 

During  that  portion  of  an  object's  trajectory  which  is  exoatmos- 
pheric,  the  object  is  subject  to  no  external  torques,  and,  therefore,  its 
motion  can  be  specified  by  the  Eulerian  angles  a,  8,  and  y  as  shown  in 
Figure  3.  For  this  to  be  true,  however,  the  angular  velocity  vector 
direction  of  the  object  had  to  have  been  chosen  as  the  z-a:.is  for  the 
body  coordinate  system.  This  axis  is  usually  the  object's  centerline. 

Movement  of  the  object  about  the  motion  axis  is  simulated  by 
describing  the  angles  a,  8,  and  y  as  functions  of  time.  This  movement 
is  referenced  to  the  motion  axis  coordinate  system  which  has  axes 
(a,  b,  MA) .  The  directions  of  the  a-  and  b-axes  will  be  specified  in 
Section  3.  To  initialize  the  angles  a,  8,  and  y  in  the  motion  axis 
coordinate  system,  assume  that  at  the  beginning  of  the  history  (t  =  0) : 

a(0)  =  aQ 

8(0)  =  6o  (4) 

y(0)  =  yo- 

The  rate  of  change  of  B(t)  determines  the  precession  of  the  object  about 
the  motion  axis  as  shown  in  Figure  3.  If  the  precession  period  is  T^, 
then  the  value  of  8  at  any  time  t  is  given  by: 

8(t)  =  s0  +  ~  .  (5) 

P 

The  value  of  y  at  the  same  time  t  is  somewhat  more  complicated. 

It  is  a  property  of  the  Eulerian  angles  that  movement  of  the  line 
of  nodes  (see  Figure  3)  in  the  a-b  plane  imparts  an  apparent  spin  on  the 
object  in  inertial  space.  More  simply,  as  the  object  with  no  spin  pre- 
cesses  about  the  motion  axis,  it  keeps  its  same  side  toward  that  axis. 

To  obtain  the  desired  inertial  spin  period  (Tg)  ,  y  must  be  calculated 
using  Equation  6. 
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FIGURE  3. 


EULER  TRANSFORMATION  FROM  THE  MOTION  AXIS  COORDINATE 
SYSTEM  TO  THE  BODY  COORDINATE  SYSTEM 
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y(t)  =  Y  +  2tt 


(Ts  '  O' 


If  nutation  is  expected,  the  angle  a(t)  Is  given  by: 


a(t)  =  ^  sin 


2TTt 


N  J 


(6) 


(7) 


where  A^  is  the  nutation  amplitude  and  is  the  nutation  period. 

The  Eulerian  angles  are  all  that  is  necessary  to  relate  the  body 
coordinate  system  to  the  motion  axis  coordinate  system.  A  vector  which 
has  a  known  direction  in  either  coordinate  system  can  be  transformed  into 
the  other  coordinate  system  using  the  D-matrix  defined  in  Equation  (A7) . 

In  Section  2.1,  it  was  shown  that  a  total  RCS  could  be  calculated  if  the 
direction  of  the  line-of-sight  to  the  radar  were  known.  If  the  LOS  vector 
is  known  in  the  motion  axis  system,  its  transformation  into  the  body  sys¬ 
tem  using  the  D-matrix  and  Equations  (5)  and  (6)  will  lead  to  a  simulation 
of  object  motion  which  can  include  spin,  precession,  and  nutation. 
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3.  TRAJECTORY  DEPENDENT  VARIABLES 


3.1  PRINCIPLE  AXIS  COORDINATE  SYSTEM 

A  number  of  objects  for  which  models  are  required  are  oriented 
relative  to  a  particular  spatial  direction  called  the  principle  axis  (PA) . 
This  axis,  which  lies  in  the  trajectory  plane,  can  be  specified  by  the 
user  or  by  flight  models  which  have  been  developed  by  TBE.  Figure  4 
shows  the  coordinate  system  which  has  been  constructed  about  the  princi¬ 
ple  axis  and  the  flight  trajectory  plane. 

The  trajectory  plane  of  an  exoatmospheric  ballistic  object  is 
formed  by  the  intersection  of  that  object's  instantaneous  velocity  vector 
and  a  vector  (gravity)  which  points  to  the  center  of  the  earth.  In  an 
inertial  frame  of  reference,  the  trajectory  plane  remains  unchanged 
during  flight  because  the  only  force  acting  on  an  exoatmospheric  ballistic 
object  is  gravity.  Since  the  trajectory  plane  is  inertially  stable,  its 
perpendicular  is  ideally  suited  for  an  inertially  fixed  coordinate  system. 

To  form  the  principle  axis  coordinate  system,  the  trajectory  plane 

A  A 

perpendicular  (_L)  and  the  principle  axis  (PA)  are  chosen  as  the  first  and 
third  axes,  respectively.  The  positive  sense  of  PA  is  such  that  the  vector 
to  the  center  of  the  earth  has  a  positive  component  along  the  Praxis,  and 
the  trajectory  plane  perpendicular  is  chosen  so  that  the  second  axis  (v) 
has  a  positive  velocity  component.  Therefore, 


The  principle  axis  coordinate  system  is  now  uniquely  defined  at  the  time 
the  principle  axis  is  chosen  and  is  invariant  during  the  exoatmospheric 
portion  of  the  flight. 

3.2  MOTION  AXIS 

The  relationship  between  the  motion  axis  and  the  principle  axis 
is  shown  in  Figure  5.  When  the  object  is  deployed  during  its  flight, 
its  motion  axis  points  in  a  particular  direction  relative  to  the  principle 
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VECTOR  TO 
CENTER  OF 
EARTH 


LOS 


1 

PERPENDICULAR  TO 
THE  TRAJECTORY  PLANE 


FIGURE  4.  THE  PRINCIPLE  AXIS  COORDINATE  SYSTEM  AND  THE 
RADAR  LINE  OF  SIGHT  (LOS) 


axis.  The  angle  between  the  P^-axis  and  the  MA-axis  is  called  PMA  (motion 
axis  polar  angle).  The  angle  between  the  trajectory  plane  perpendicular 
and  the  projection  of  the  motion  axis  in  the  j_-v  plane  is  called  AMA 
(motion  axis  azimuth  angle) . 

A  new  coordinate  system  —  the  motion  axis  coordinate  system  — 
is  now  defined  with  MA  as  its  third  axis  and  a  and  b  as  its  first  and 
second  axes,  respectively,  as  shown  in  Figure  5.  The  transformation 
matrix  from  the  principle  axis  coordinate  system  to  the  motion  axis 
coordinate  system  is  defined  by  the  A-matrix  described  in  the  appendix. 

The  directions  a  and  b  are  uniquely  defined  by  the  transformation  in 

A 

Equation  (A5) .  If  the  object  is  deployed  so  that  the  MA-axis  aligns 
exactly  with  the  principle  axis,  the  A-matrix  becomes  the  identity  matrix. 

A 

It  was  stated  in  Section  2.2  that  if  the  LOS-vector  were  known 
in  the  motion  axis  coordinate  system,  an  RCS  calculation  could  be  made. 

/s 

Continuing  the  logic,  if  the  LOS  vector  is  known  in  the  principle  axis 
coordinate  system,  it  can  be  transformed  into  the  motion  axis  system 
using  the  A-matrix  defined  by  Equation  (A5) .  The  characterization  of 

A 

the  LOS  vector  in  the  principle  axis  coordinate  system  is  discussed  in 
the  following  section. 


A.  VIEWING  DEPENDENT  VARIABLES 


As  the  object  moves  along  its  flight  path,  the  direction  of  the 
vector  from  the  center-of-gravity  (origin  of  the  principle  axis  coordinate 

A 

system)  to  the  radar  location  (LOS)  changes  —  see  Figure  4.  Thus,  in  the 
principle  axis  system,  0 ( t )  and  0(t)  are  functions  of  time.  The  time 
history  of  these  two  angles  defines  the  inertial  direction  of  the  principle 
axis,  characterizes  the  flight  trajectory  and  the  trajectory  plane,  and 
pinpoints  the  radar  location.  The  TBE  developed  flight  models  can  be 
used  to  generate  0(t)  and  <f)(t)  for  any  trajectory,  principle  axis,  and 
radar  location  for  use  in  the  MSPB  motion  generator.  These  calculations 
can  be  made  in  real  time  during  a  simulation,  but  time  can  be  saved  using 
other  techniques. 

If  the  viewing  parameters  are  such  that  the  angle  functions  6(t) 
and  <}>(t)  are  reasonably  smooth,  these  angles  can  be  approximated  by  an 
Nth  power  polynominal  in  time.  Essentially: 

N0 

0(t)  =  l  An  tn  (8) 

n=0 


<J>(t)  =  l  B  t  • 
n=0 


(9) 


Should  they  not  be  smooth,  the  angles  can  be  listed  in  a  look-up  table  as 
a  function  of  time.  In  either  case,  a  real  time  simulation  can  be 
implemented . 

/\ 

Once  9(t)  and  $(t)  are  known  at  a  given  moment,  the  LOS-vector 
in  the  principle  axis  coordinate  system  is  given  by: 

LOS  (j_,  v,  P)  =  _L  sin  9  cos  <j>  +  v  sin  9  sin  <p  +  J_  cos  0  .  (10) 

A 
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5.  SUMMARY 


In  essence,  all  motion  is  relative  to  the  principle  axis 
coordinate  system  which  does  not  move.  As  the  object,  properly  oriented, 
flies  along  its  chosen  trajectory,  the  apparent  motion  of  the  viewing 
radar  is  modeled  b>  Equations  (8),  (9),  and  (10).  The  viewing  dependent 
parameters  are  listed  in  Table  1.  The  time  history  of  the  line-of-sight 
motion  in  the  principle  axis  system  is  the  first  step  in  the  modeling 
which  is  aspect  dependent. 

Since  the  motion  axis  coordinate  system  is  not  generally  aligned 
with  the  principle  axis  coordinate  system  —  namely  the  angular  ,iomo  ltum 
vector  of  the  object  is  not  colinear  with  the  principle  axis  —  c  tans- 
formation  between  the  two  coordinate  systems  is  necessary.  That 

A 

transformation  is  Equation  (A5)  anti  is  used  to  transform  the  LOS-vector 
into  the  motion  axis  system. 

LOS  (a,  b,  MA)  =  A  •  LOS  (I,  v,  Pj  (11) 

The  parameters,  called  trajectory  dependent  variables,  associated  with 
this  transformation  are  listed  in  Table  2.  The  range  of  values  of  these 
variables  depends  on  the  object  being  modeled  and  somewhat  on  the 
trajectory  along  which  the  object  is  moving. 

Object  spin,  precession,  and  nutation  are  simulated  in  the 
transformation  from  the  motion  axis  .  ordinate  system  to  the  body 
coordinate  system  as  defined  in  Equation  (A7) .  These  parameters,  which 
are  listed  in  Table  3,  are  essentially  functions  only  of  the  object 
being  modeled,  and  their  ranges  cannot  be  set  without  knowing  something 
of  that  object.  Additional  parameters  may  become  necessary  which  will 
complete  the  model.  The  LOS-vector  in  the  body  system  is  then 

LOS  (x,  y,  z)  =  D  •  LOS  (a,  b,  MA)  .  (12) 

Substituting  Equation  (11)  into  (12)  yields  the  final  equation  for 
modeling  the  motion  of  an  object  flying  along  its  trajectory. 
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TABLE  1.  MSPB  VIEWING  DEPENDENT  VARIABLES* 


SYMBOL 

DESCRIPTION 

An  n  =  0,  1,  3  . 

■ 

Coefficients  of  the  polynominals 
describing  the  movement  of  the 

B_  n  =  0,  1,  3  . 

- 

1 ine-of-sight  in  the  principle 

1 

axis  coordinate  system 

The  angles  e(t)  and  could  also  be  calculated  directly  from 
a  trajectory. 


TABLE  2.  MSPB  TRAJECTORY  DEPENDENT  VARIABLES 


SYMBOL 

DESCRIPTION 

A 

P 

A 

The  direction  of  the  principle 
axis 

PMA 

Angle  between  the  principle 
axis  and  the  motion  axis 

(radians) 

AMA 

Angle  between  the  trajectory 

plane  perpendicular  and  the 
projection  of  the  motion  axis 
into  the  ±-v  plane  (radians) 

TABLE  3.  MSPB  TRAJECTORY  INDEPENDENT  VARIABLES 


SYMBOL 

DESCRIPTION 

<V  60’  Yo 

Initial  coning,  precession,  and 
spin  angles  (radians) 

Tp 

Precession  period  (seconds) 

Ts 

Spin  period  (seconcs) 

tn 

Nutation  period  (seconds) 

an 

Nutation  amplitude  (radians) 
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6.  EXAMPLE  IMPLEMENTATION 


MSPB  motion  parameters  have  been  defined  in  the  preceeding  sections. 
These  are  listed  in  Tables  1,  2,  and  3.  The  use  of  the  MSPB  generator 
requires  a  definition  of  each  parameter.  In  general,  these  definitions 
will  be  in  the  form  of  probability  density  functions  to  account  for  the 
variability  in  object  motion  that  is  typically  observed  on  various  test 
flights  of  the  same  object. 


To  complete  the  MSPB  model,  the  PCS  (0^)  and  position  (S^)  of 
each  scattering  components  is  required.  As  an  example,  consider  an 
object  with  two  components: 


CJj  =  16 


sin  x 


x  =  23  sin  6^ 


(14) 


=  16 


I  sin  Y 


Y  =  8  cos  <J> 

L 


(15) 


Sx  =  £  x  (18) 

V=  0  .  (18) 

A 

Scatterer  #1  has  an  RCS  which  is  symmetric  about  the  z-axis  and  is  not 
located  at  the  center  of  motion  but  rather  some  distance  along  the 
positive  x-axis.  Scatterer  #2  is  located  at  the  center  of  motion  but  is 
not  symmetric  about  the  z-axis.  The  RCS  static  patterns  of  these  two 
scattering  components  are  shown  in  Figures  6  and  7.  Assume  that  the 
phase  delay  (A^)  is  zero  for  each  component. 

Figures  8  through  11  show  examples  of  the  MSPB  time  varying  RCS 
output  using  the  scatterers  described  above  in  40-second  histories 
sampled  at  a  24  Hz  data  rate.  The  MSPB  parameters  are  shown  in  Table  4. 


Figure  8  shows  the  effect  of  precession  only  while  Figure  9  is  the  result 
of  spin  only.  Both  figures  show  the  effect  of  moving  along  the  flight 
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29. 
1  S. 
1  9. 
5. 
9. 
-5. 
-to. 
-I  5. 
-29. 
-25. 
-39. 
-35. 
-49. 


FIGURE  6.  RCS  STATIC  PATTERN  OF  SCATTERER  #1 
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I  5. 
1  0. 


S. 


0. 


-S. 


-1  0. 


-IS. 

-20. 


-25. 


FIGURE  8.  EFFECT  OF  PRECESSION  ONLY  ON  THE  EXAMPLE  OBJECT 


TIME  (seconds) 


FIGURE  9.  EFFECT  OF  SPIN  ONLY  ON  THE  EXAMPLE  OBJECT 


TABLE  4.  MSPB  PARAMETERS  USED  IN  THE  EXAMPLE  IMPLEMENTATION 


SYMBOL 

VALUE 

A,. 

2.0 

A2 

0.006 

A3 

0.00005 

Bi 

15.0 

b2 

-0.7 

8  3 

0.003 

A 

PA 

Vertical 

PMA 

0.07  (4  degrees) 

AMA 

*/2 

a 

0 

0.05  (3  degrees) 

8o’  Y 0 

0 

TP 

10  sec 

Ts 

2  sec 

an 

0.07  (4  degrees) 

tn 

2  sec 
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path  which  has  been  modeled  as  the  2n<^  order  polynominal  in  time.  The 
combined  effects  of  both  spin  and  precession  is  shown  in  Figure  10,  and 
a  small  nutation  has  been  added  in  the  data  shown  in  Figure  11. 

Thus  this  model  is  completely  specified  by  the  trajectory  depen¬ 
dent  and  independent  parameters  of  Tables  2  and  3,  the  main  body  RCS 
static  pattern  and  an  algorithm  for  the  scatterer  RCS.  If  the  scatterer 
is  allowed  to  move  relative  to  the  main  body,  a  model  of  this  relative 
motion  is  also  required.  The  method  for  determining  the  composite  RCS 
must  also  be  specified  for  any  given  model. 
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APPENDIX  A.  DERIVATION  OF  THE  TRANSFORMATION  MATRICES  (A  AND  D) 


If  two  coordinate  systems  with  a  common  origin  are  defined  in 
space  such  that  their  respective  axes  do  not  coincide,  the  problem  often 
arises  of  calculating  the  components  of  a  vector  in  one  system  from  the 
vector  components  as  described  in  the  second  system.  Such  a  conversion 
of  vector  components  is  called  a  transformation  and  is  best  described 
by: 


!'  =  l  a. .  V. 

1  j=l  ^  J 


i  -  1,  2,  3. 


(Al) 


The  quantities  (Vj,  V2,  V3)  are  defined  as  the  components  of  the  vector 

— ►  — ►  ^ 

V  in  the  unprimed  coordinate  system.  The  second  vector  V  is  defined  in 

terms  of  the  transformation  properties  of  its  components  under  a  rotation’ 
of  the  coordinate  system.  The  a_  are  identified  as  the  direction  cosines 
of  the  angle  as  measured  from  the  i-axis  in  the  primed  system  to  the 
j-axis  in  the  unprimed  system.  A  rectangular  coordinate  system  can 
always  be  brought  into  coincidence  with  a  second  rectangular  coordinate 
system  with  the  same  origin  by  rotating  (transforming)  about  each 
orthogonal  axis  in  turn  and  then  combining  the  results.  This  informa¬ 
tion  will  now  be  used  to  derive  the  transformation  matrices  (A  and  D) 
used  in  Sections  3  through  5. 

The  transformation  from  the  principle  axis  system  to  the  motion 
axis  system  involves  first  a  rotation  of  AMA  degrees  about  the  P-axis 
(see  Figure  5).  The  direction  cosines  for  that  rotation  are: 

bn  =  cos  (90  -  AMA)  =  sin  (AMA) 
bi2  *  cos  (180  -  AMA)  -  -cos  (AMA) 
bn  =  cos  (90)  =  0 
b2i  =  cos  (AMA) 

b 22  =  cos  (90  -  AMA)  -  sin  (AMA)  (A2) 

b 23  =  cos  (90)  =  0 
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bai 

=  COS 

(90)  =  0 

b  32 

=  cos 

(90)  =  0 

b  33 

*  cos 

(0)  =  1. 

A  second  rotation  of  PMA  degrees  about  a-axis  will  bring  the  principle 
axis  system  in  line  with  the  motion  axis  system. 


Cll 

= 

COS 

(0)  = 

1 

C  12 

= 

COS 

(90) 

=  0 

C13 

= 

COS 

(90) 

=  0 

C21 

= 

COS 

(90) 

=  0 

C22 

= 

COS 

(PMA) 

C23 

= 

COS 

(90  +  PMA)  = 

C  31 

= 

COS 

(90) 

=  0 

C  32 

S 

COS 

(90  - 

PMA)  = 

C33 

= 

COS 

(PMA) 

# 

=  -sin  (PMA) 


=  sin  (PMA) 


(A3) 


The  total  transformation  is  given  by  multiplying  Equations  (A2)  by  (A3) . 
Namely: 


a±j  k=!  Cik  bkJ 


i  =  1,  2,  3 


1,  2,  3. 


(A4) 


The  aJ ,  are  the  elements  of  transformation  matrix  A  which  is  then: 
ij 

(sin  (AMA)  -cos  (AMA)  0 

cos  (PMA)  cos  (AMA)  cos  (PMA)  sin  (AMA)  -sin  (PMA) j(A5) 

sin  (PMA)  cos  (AMA)  sin  (PMA)  sin  (AMA)  cos  (PMA)/ 


A 


Any  vector  V  can  then  be  transformed  from  the  principle  axis  coordinate 
system  (unprimed)  to  the  motion  axis  coordinate  system  (primed)  using: 


V'  =  A  •  V 

which  is  the  equivalent  of  Equation  (Al) . 


(A6) 
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The  transformation  shown  in  Figure  3  is  a  standard  Euler  trans¬ 
formation*.  The  D-transformation  is  a  rotation  first  about  the  MA-axis, 
second  about  the  line  of  nodes,  and  finally  about  the  new  z-axis.  The 

elements  of  the  Q-matrix  which  will  carry  out  the  transformation  from 

the  motion  axis  system  to  the  body  coordinate  system  are  listed  in 
Equation  (A7) . 

dn  =  cos  y  *  cos  8  -  cos  a  •  sin  8  *  sin  y 

di2  =  cos  y  *  sin  8  +  cos  a  •  cos  8  *  sin  y 

dn  =  sin  y  •  sin  a 

d2i  =  -sin  Y  *  cos  8  -  cos  a  •  sin  8  *  cos  y  (A7) 

d22  =  cos  a  •  cos  8  •  cos  y  -  sin  y  *  sin  8 

d23  =  cos  y  •  sin  a 

d3i  =  sin  a  •  sin  8 

d32  =  -sin  a  •  cos  8 

d 33  =  cos  a. 


In  the  body  coordinate  system  (double-primed)  the  V'  is  then: 


V 


i 


3 

I 


k=l 


dik  Vk 


(A8) 


or 

V"  =  D  •  V".  (A9) 

The  value  of  V'  [Equation  (A6) ]  can  be  substituted  into  Equation  (A9)  to 
yield  a  direct  transformation  from  the  principle  axis  coordinate  system 
to  the  body  coordinate  system. 

V"  =  e  •  4  •  V  .  (A10) 


*Goldstein,  H.,  "Classical  Mechanics",  Addison-Wesley ,  1965,  pp.  107-109. 
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The  vector  V  has  components  which  are  properly  defined  in  the 
(x,  y,  z)  coordinate  system  shown  in  Figure  2.  Namely: 


V"  »Vx  +  Vy  +  Vz. 
x  y  z 


(All) 


Therefore  the  polar  angle  (9  )  and  the  azimuth  angle  (<b  )  are  then 

L  L 

given  by: 


0L  «  cos-1  (Vz) 


and 


(A12) 


(A13) 
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APPENDIX  C 


CORRELATED  STATISTICAL  (CORS)  TIME  VARYING  RCS  ALGORITHM 


By 

W.  J.  Gruner 
L.  K.  Montgomery,  Jr. 


C-l 


1.  INTRODUCTION 


TBE  is  required  to  develop  exoatmospheric  time  varying  L-band 
and  UHF  RCS  models  for  objects  of  several  weapons  systems  under  Contract 
DASG60-79-C-0094.  A  statistical  generation  method  is  one  procedure  for 
simulating  RCS  histories  that  do  not  have  significant  time  structure. 

TBE  has  developed  a  correlated  statistical  time  varying  RCS  model  for 
this  purpose. 

Section  2  presents  a  description  of  the  correlated  statistical 
time  varying  RCS  model.  This  model  has  been  used  with  Gaussian  probability 
density  functions  and  the  parameters  for  this  form  of  the  model  are  given. 
Representative  time  histories  are  presented  for  several  sets  of  input  para¬ 
meters. 
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2.  CORRELATED  STATISTICAL  RCS  GENERATION  MODEL 


The  correlated  statistical  RCS  generation  model  provides  an  RCS 
value  given  the  last  RCS  value  and  the  intervening  time  interval  (At) . 

It  is  represented  mathematically  by 

RCS  =  p*RCS  .  ,  +  (1-p)  R 
new  old 

where  p  =  exp  (-COR  *At)  =  correlation  coefficient 

and  R  =  new  random  draw. 

The  exponential  correlation  function  was  based  on  the  typical  exponential 
decay  observed  in  autocorrelation  functions  of  RCS  time  histories.  The 
underlying  probability  density  is  based  on  the  observed  data  character¬ 
istics.  Without  additional  constraints,  this  procedure  can  give  unrealis¬ 
tic  RCS  values  since  mathematically  the  limits  of  a  probability  density 
may  extend  beyond  physical  reality.  Thus  the  model  provides  a  limit  to 
the  amount  the  RCS  can  change  in  a  time  interval  and  both  upper  and  lower 
RCS  limits. 

A  flow  diagram  of  this  model  is  shown  in  Figure  1.  The  random 
number  R  is  drawn  from  the  specified  probability  density.  The  bound  is 
the  maximum  change  in  RCS  allowed  in  the  time  interval  At.  ant* 

RCS>hn  are  specified  to  keep  the  RCS  values  within  the  realm  of  physical 
reality.  When  these  limits  are  exceeded,  the  new  RCS  value  is  changed 
to  stay  within  the  specified  limits. 

The  initial  use  of  this  model  has  been  with  the  probability  density 


defined  as  a  Gaussian  distribution  in  dBsm  space.  For  the  correlated 

Gaussian  time  varying  RCS  model,  the  six  required  parameters  are  given 

in  Table  1.  It  was  found  to  be  convenient  to  define  BOUND,  RCS.,.„  and 

_ _ MAX _ 


SYMBOL 

parameter  definition 

U 

Mean  of  Gaussian  Density 

a 

Standard  Deviation  of  Gaussian  Density 

COR 

Correlation  Coefficient 

BOUND 

Maximum  RCS  Change  in  At 

rcsmax 

Maximum  Allowed  RCS  Value 

rcsmin 

Minimum  Allowed  RCS  Value 

TABLE  1.  CORRELATED  GAUSSIAN  TIME  VARYING  RCS  MODEL  PARAMETERS 
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FIGURE  1.  CORRELATED  STATISTICAL  TIME  VARYING  RCS  MODEL 
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RCSiitv(  in  terms  of  p  and  0  as 
MIN 

BOUND  =  NUM*a*exp  (At) 

RCSMAX  “  M  +  NSU*° 
rcsmin  “  w  -  NSL*°- 

Thus  for  the  correlated  Gaussian  model,  the  six  parameters  are  p,  0, 

COR,  NUM,  NSU,  and  NSL.  To  model  a  particular  object,  these  parameters 
will  in  general  be  defined  as  probability  distributions. 

In  developing  this  model  it  was  found  necessary  to  expand  the 
Gaussian  density  for  highly  correlated  draws  (i.e.  small  values  of  COR 
and  At).  This  expansion  is  accomplished  by  using  a1  =  0/  FACTOR  instead 
of  O  as  the  standard  deviation  of  the  Gaussian  distribution  where 

/  1  ;  C0R*At S  2.86 

FACTOR  =  l  0.6056*logi()  (C0R*At)  +  0.723;  otherwise 

Examples  of  the  model  output  for  selected  values  of  the  parameters 
are  shown  in  Figure  2.  These  data  were  generated  at  24  Hz.  More  struc¬ 
ture  is  seen  in  Figure  2(c)  which  has  the  smallest  correlation  coefficient 
and  standard  deviation  of  the  cases  illustrated. 
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APPENDIX  D 


SINUSOIDS  AND  RANDOM  PROCESS  (SARP) 
TIME  VARYING  RCS  ALGORITHM 


By 

W.  J.  Gruner 


1.  INTRODUCTION 


TBE  is  required  to  develop  exoatmospheric  time  varying  L-band 
and  UHF  RCS  models  for  objects  of  several  missile  systems  under  Contract 
DASG60-79-C-0094 .  TBE  has  developed  SARP  for  the  simulation  of  RCS  time 
histories  from  objects  with  significant  time  structure  and  without 
definable  motion  characteristics.  SARP  provides  simulated  RCS  values  by 
the  sum  of  harmonically  related  sinusoids  and  an  uncorrelated  random 
process  in  dBsm  space. 

Section  2  presents  a  description  of  SARP.  Model  parameters, 
generation  variables,  and  constraints  are  discussed.  Representative  time 
histories  are  presented  for  several  sets  of  input  parameters. 
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2.  RCS  GENERATION  MODEL  SARP 


SARP  provides  an  RCS  value  at  a  specified  time  given  sinusoidal 
amplitudes,  frequencies,  initial  phases,  and  relative  random  process 
components.  A  diagram  of  SARP  is  shown  in  Figure  1.  The  returned  RCS 
value  in  dBsm  is  expressed  mathematically  as: 

RCS  —  p  +  Rg  +  Rjj 

where 

U  is  the  desired  mean  value 

Rs  is  the  sinusoidal  component 
'  from  NSIN  harmonically  related  sinusoids, 

Rjj  is  a  Gaussian-distributed  random  component.* 

These  variables  are  functionally  related  to  the  required  SARP  model 
parameters  in  Table  1.  Some  of  these  parameters  -  M,  VAR,  RELN,  and  F0  - 
would  generally  be  specified  in  terms  of  probability  density  functions 
to  provide  the  variability  between  various  observations  of  the  same  object. 
The  number  of  sinusoids  (NSIN),  relative  amplitudes  (AMPN)  and  initial 
phase  (PHAS)  are  determined  from  power  spectral  density  analyses  of 
measured  RCS  history  data. 

The  SARP  random  variable  Rjj  is  a  zero  mean  Gaussian  random 
variable*  with  standard  deviation  aN.  Since  mathematically,  the  limits  of 
this  distribution  can  extend  beyond  physical  reality,  the  value  of  R^  is 
constrained  by  |r^|  _<  30jj.  The  standard  deviation  is  related  to  the 
parameters  RELN  and  VAR  by 

„„  -  10  WLN/10)  ,  VAR  >S 

where  R^j  and  0^  have  the  units  of  dBsm.  The  relative  random  process 
component  (RELN)  is  constrained  by  RELN  <_  0.  If  RELN  =  0  the-e  is  no 

*  If  appropriate,  other  probability  density  functions  could  be  used  for  Rjj. 
However,  all  applications  to  date  by  TBE  have  used  Gaussian  probability 
density  functions  so  the  model  is  described  in  this  manner. 


D-3 


1 


FIGURE  1.  TIME  VARYING  RCS  GENERATION  MODEL  SARP 
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TABLE  1 .  SARP  TIME  VARYING  RCS  MODEL  PARAMETERS 


SYMBOL 

PARAMETER 

u 

Process  mean  (dBsm) 

VAR 

Process  variance  (dBsm2) 

RELN 

j 

Relative  random  process  component 
(dB  relative  to  composite) 

F_ 

Sinusoidal  component  fundamental  frequency 

0 

(Hz) 

NSIN 

Number  of  harmonically  related  sinusoids 
(unitless) 

AMPN 

Array  of  sinusoidal  amplitudes  normalized 
relative  to  sinusoidal  composite 
(unitless) 

PHAS 

Array  of  sinusoidal  initial  phase  (deg) 

TABLE  2.  EXAMPLE  SARP  SINUSOIDAL  INPUT,  NSIN  = 


.10 


-140 


sinusoidal  component  and  the  simulation  consists  of  a  constrained 
correlated  random  process. 


un- 


The  SARP  sinusoidal  component  R  is  determined  bv 

s  J 

NSIN 

R  =  l  AMP  *  SIN  (FREQ,  *  (T  -T  ,,)  +  PHAZJ 
*  ±  =  l  x  1  otr  i 

Where  AMPi  +  AMPNi  *  (VAR  *  1.  _10RELN/10  ) 

is  the  sinusoidal  amplitude  (dBsm) 

FREQ,  =  i  *  F  *  2tt 
1  o 

is  the  sinusoidal  frequency  (rad/sec) 

To££  is  a  uniform  initial  bias  time  drawn  from  the 

interval  (0,  1/F  ) 
o 

PHAZ  =  (PHAS) 

is  the  sinusoidal  initial  phase  (rad)  at  time  T  =  0. 

Physical  constraints  are  modeled  by  limiting  the  amplitude  densities 
associated  with  u,  VAR,  RELN  and  FQ. 

Examples  of  SARP  output  for  selected  input  parameters  are  shown 
in  Figure  2.  These  figures  were  generated  with  NSIN  =  4  and  values  for 
relative  amplitudes  and  initial  phase  angle  are  given  in  Table  2. 

These  data  were  generated  at  the  indicated  data  rates.  More 

structure  is  seen  as  RELN  is  decreased.  Figures  2(a)  -  (c)  are  40-second 

realizations  sampled  at  a  constant  data  rate  of  24  Hz.  Figure  2(a)  is 

predominantly  sinusoidal  with  RELN  =  -20dB  and  F0  =  3.31  Hz.  Figures  2(b) 

and  2(c)  were  generated  with  F0  =  1.83  Hz  and  RELN  =  -10  and  -2.  With 

RELN  =  -2  the  structure  of  the  simulated  pattern  is  predominantly  random. 

Figures  2(d)  -  (f)  are  200-second  realizations  of  the  same  processes  as 

used  in  Figures  2(a)  -  (c)  but  with  much  lower  data  rates  of  1  Hz  in 

Figure  2(d)  and  0.5  Hz  in  Figures  2(e)  and  2(f).  In  these  figures  the 

sinusoidal  pattern  is  also  shifted  by  T  ,  =  0.5  sec  at  simulation  time  zero. 

or  t 

Figure  3  shows  SARP  realizations  at  time  intervals  corresponding 
to  non-uniform  data  rates  obtained  on  an  object  for  which  SARP  has  been 
used.  The  parameters  incorporated  in  SARP  for  Figures  3(a)  and  3(b)  are 


D-6 


identical.  With  the  relative  random  noise  (RELN)  at  -4  dB,  the  sinusoidal 
and  the  random  components  contribute  approximately  60%  and  40%  in  terms  of 
RCS  variance.  The  average  data  rate  in  Figure  3(a)  is  approximately  15  Hz 
over  an  interval  of  nearly  40  seconds  and  the  average  data  rate  in  Figure 
3(b)  is  0.67  Hz  over  an  interval  of  approximately  170  seconds.  In  the 
figures,  the  irregular  sampling  causes  the  simulated  histories  to  appear 
predominantly  random.  In  Figure  3(b),  if  a  sinusoidal  process  were  postu¬ 
lated,  an  observer  might  deduct  an  apparent  frequency  of  0.1  Hz  as  opposed 
to  the  true  SARP  fundamental  periodicity  associated  with  3.31  Hz. 
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FIGURE  2.  EXAMPLE  SARP  RCS  REALIZATIONS,  SAMPLING  AT  CONSTANT  DATA  RATES 
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1.  INTRODUCTION 


TBE  is  required  to  develop  exoatmospheric  time  varying  L-band 
and  UHF  RCS  models  for  objects  of  several  missile  systems  under  Contract 
DASG60-79-C-0094.  TBE  has  developed  MSARP  for  the  simulation  of  RCS  time 
histories  from  complex  objects  with  significant  time  structure  and  with¬ 
out  definable  motion  characteristics.  MSARP  provides  simulated  RCS  values 
by  the  sum  of  related  sinusoids  and  an  uncorrelated  random  process  in 
dBsm  space  as  provided  by  SARP  (Appendix  D) .  MSARP  expands  the  concept  of 
SARP  by  randomly  selecting  higher  frequency  components. 

Section  2  presents  a  description  of  MSARP.  Model  parameters, 
generation  variables,  and  constraints  are  discussed.  Representative  time 
histories  are  presented  for  several  sets  of  input  parameters. 
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2,  RCS  GENERATION  MODEL  MSARP 

MSARP  provides  an  RCS  value  at  a  specified  time  based  on  sinusoidal 
amplitudes,  frequencies,  initial  phases,  and  relative  random  process 
components.  A  diagram  of  MSARP  is  shown  in  Figure  1.  The  returned  RCS 
value  in  dBstn  is  expressed  mathematically  as: 

RCS  =  p  +  R„  +  R 
“  N 

where 

U  is  the  desired  mean  value, 

Rg  is  the  sinusoidal  component 

from  NSIN  sinusoids, 

R^  is  a  Gaussian-distributed  random  component. 

These  variables  are  functionally  related  to  the  required  MSARP  model 
parameters  in  Table  1.  Some  of  these  parameters  (p,  VAR,  RELN,  and  Fq) 
would  generally  be  specified  in  terms  of  probability  density  functions 
to  provide  the  variability  between  various  observations  of  the  same  object. 

The  MSARP  random  variable  R^  is  a  zero  mean  Gaussian  random  vari¬ 
able*  with  standard  deviation  oN-  Since  mathematically,  the  limits  of 
this  distribution  can  extend  beyond  physical  reality,  the  value  of  R^  is 
constrained  by  |Rjj|<3aN.  The  standard  deviation  is  related  to  the 
parameters  RELN,  VAR,  and  Rg  by 

O  =  MAX  JO. 01, 1.8  (10  (RELN/1°)  *  varS  (0.85  -  0.04  *  R  ) ! 

where  R^  and  0^  have  the  units  of  dBsm.*  The  relative  random  process 
component  (RELN)  is  constrained  by  RELN  s  0.  If  RELN  =  0  there  is  no 
sinusoidal  component  and  the  simulation  consists  of  a  constrained  uncorre¬ 
lated  random  process. 

For  MSARP,  the  number  of  sinusoids  is  randomly  chosen  from  a  uni¬ 
form  density  with  a  specified  range.  This  range  is  derived  from  power 
spectral  density  analyses  of  the  associated  data.  The  sinusoidal  frequen¬ 
cies  are  chosen  to  be  nonharmonic  by  a  small  random  component.  For  example. 


*If  appropriate,  other  probability  densisy  functions  could  be  used  for  R, 
However,  all  applications  to  date  by  TBE  have  used  Gaussian  probability 
density  functions  so  the  model  is  described  in  this  manner. 


FIGURE  1.  TIME  VARYING  RCS  GENERATION  MODEL  MSARP 
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TABLE  1.  MSARP  TIME  VARYING  RCS  MODEL  PARAMETERS 


SYMBOL 

PARAMETER 

P 

Process  mean  (dBsm) 

VAR 

Process  variance  (dBsm2) 

RELN 

Relative  random  process  component 

(dB  relative  to  composite) 

F 

Sinusoidal  component  fundamental  frequency 

0 

(Hz) 

NSIN 

Number  of  sinusoids  (unitless) 

AMPN 

Array  of  relative  sinusoidal  amplitudes 

for  initial  subset  of  sinusoids  (unitless) 

0 

Standard  deviation  of  frequency  non- 

F 

harmonic  component 

a 

Standard  deviation  of  sinusoidal  amplitude 

A 

random  component 

Logical  experession  for  amplitudes  of  sinusoids  of  higher  order  than 

initial  subset,  relative  to  final  specified  amplitude.  This  expres- 

sion  is  prior  to  application  of  random  component.  For  example 

AMPN ( I )  =  AMPN(J)  (1. 0-1.1  log1Q  ( I/J ) ) ,  J<I<_NSIN 

where  J  sinusoids 

were  specified  in  the  initial  subset 

where 


F.  =  f  *  (l  +  Fj  *  i 

1  O  \  til 

F..  is  a  zero-mean  Gaussian  distributed  random  variable 
N 

with  a  specified  standard  deviation,  a  ,  constrained 

r 

to  be  within  ±3o_. 

F 

Relative  sinusoidal  amplitudes  for  an  initial  subset  of  sinusoids  are 
determined  directly  from  the  power  spectral  density  analyses.  The  power 
spectral  densities  are  also  used  to  derive  a  logical  relationship  for 
establishing  the  amplitude  of  the  sinusoids.  This  logical  relationship 
is  determined  in  the  form  of  average  relative  amplitude  vs  frequency. 

The  actual  amplitude  is  then  derived  as  the  deterministic  average  ampli¬ 
tude  plus  a  random  component  where  amplitude  is  also  determined  from  power 
spectral  density  analysis.  For  example,  if  8  sinusoidal  amplitudes  have 
been  determined  the  amplitudes  of  the  remaining  NSIN  -  8  sinusoids  may  be 
determined  by 

AMPN.  =  MAX  0-01,  AMPA*AMPNg  x  (1.0-1. 1  log1()(i/8))' 

,  A  MO  A  MAvf1-  +  AMPSI6I 

where  AMPA  =  MAX  , 

(  U‘  ) 

and  AMPSI6  is  a  zero  mean  Gaussian  random  variable  with  standard  deviation  a 
with  amplitude  constrained  to  3  After  all  relative  amplitudes  have 

been  calculated,  they  are  normalized  to  produce  the  desired  variance.  For 
MSARP ,  all  initial  phase  angles  (PHA2.)  are  chosen  from  a  uniform  densitv 
ranging  from  0  to  2rr  radians. 

The  MSARP  sinusoidal  component  Rg  is  determined  identically  to  the 
SARP  technique  by 

NSIN 

R  =  Z  AMP  *  SIN  (FREQ.  *  (T-T  . J  +  PHAZ  ) 
s  ^  ^  i  i  off  iy 

where  AMP^  =  AMPN.  *  (VAR  *  jl.  -10RELN^10} )'5 

is  the  sinusoidal  amplitude  (dBsro) 

f  ~  ^ s  a  uniform  initial  time  bias  drawn 
from  the  interval  (0,  1/FQ) 

F.-6 


FREQ.  =  F.  *  2tt 
i  l 

is  the  sinusoidal  frequency  (rad/sec) 

Physical  constraints  are  modeled  by  limiting  the  amplitude  densities 
associated  with  y,  VAR,  RELN  and  F  . 

Example  HSARP  outputs  for  the  parameters  of  Table  2  are  shown  in 
Figure  2.  All  realizations  were  derived  at  constant  values  of  y,  VAR, 
and  RELN.  Figures  2  (a-c)  were  generated  at  a  fundamental  frequency  of 
0.025  Hz  and  Figures  2  (d-f)  at  a  fundamental  frequency  of  0.010  Hz.  All 
data  are  formulated  for  a  period  of  80  seconds  at  a  constant  20  Hz  data  rate 
These  examples  are  repetative  except  for  the  random  component  in  the 
period  associated  with  the  fundamental  frequency,  and  can  be  used  to 
simulate  complex  scatters. 
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TABLE  2.  EXAMPLE  MSARP  INPUT 
(INITIAL  SINUSOIDAL  INPUT  FOR  8  FREQUENCIES) 


SINUSOIDAL 

AMPN 

1 

2 

3 

n 

5 

6 

a 

8 

.76 

.53 

.41 

.45 

.42 

.54 

.23 

.50 

AMPN(I)  =  AMPN  (8)  •  (1.0  -  1.1  1og1Q 

(1/8)),  I.  >  8 

a 

A  =  0.2 

a  =  0.0016 

F 
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(d)  20  Hz  Data  Rate 


(e)  20  Hz  Data  Rate 
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